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Abstract 

Motivated by the dark energy issue, the one-loop quantization approach for a family of 
relativistic cosmological theories is discussed in some detail. Specifically, general f{R) gravity 
at the one-loop level in a de Sitter universe is investigated, extending a similar program 
developed for the case of pure Einstein gravity. Using generalized zeta regularization, the 
one-loop effective action is explicitly obtained off-shell, what allows to study in detail the 
possibility of (de)stabilization of the de Sitter background by quantum effects. The one-loop 
effective action maybe useful also for the study of constant curvature black hole nuclcation 
rate and it provides the plausible way of resolving the cosmological constant problem. 
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1 Introduction 



Recent astrophysical data indicate that our universe is currently in a phase of accelerated ex- 
pansion. One possible explanation for this fact is to postulate that gravity is being nowadays 
modified by some terms which grow when curvature decreases. This could be, for instance, an 
inverse curvature term ^ which might have its origin at a very fundamental level, as string/M- 
theoryjjj or the presence of higher dimensions [3^, which seem able to explain such accelerated 
expansion. Gravity modified with inverse curvature terms is known to contain some instabilities 
^ and cannot pass some solar system tests, but further modifications of the same which include 
higher derivative, curvature squared terms make it again viable (The Palatini formulation 
may also improve the situation, for recent discussion, see and refs. therein). 

Having in mind possible applications of modified gravity for the late time universe, the 
following question appears. If it so happens that Einstein gravity is only an approximate the- 
ory, looking at the early universe should this (effective) quantum gravity be different from the 
Einsteinian one? A widely discussed possibility in this direction is quantum gravity (for 
a review, see [Tj). However, other modifications are welcome as well, because they sometimes 
produce extra terms which may help to realize the early time inflation. This is supported by 
the possibility of accelerated expansion with simple modified gravity. Thus, we will study here 
general f{R) gravity at the one-loop level in a de Sitter universe. A similar program for the 
case of pure Einstein gravity (recall that it is also multiplicatively non-renormalizable) has been 
initiated in refs. |H1 IHl (see also JJ). Using generalized zeta- functions regularization (see, 
for instance |12| I13j). one can get the one-loop effective action and then study the possibility 
of stabilization of the de Sitter background by quantum effects. Moreover, such approach hints 
also to a possible way of resolving the cosmological constant problem jlO,,. Hence, the study of 
one- loop f{R) gravity is a natural step to be undertaken for the completion of such a program, 
keeping always in mind, however, that a consistent quantum gravity theory is not available yet. 
But in any case, one should also not forget that, from our present knowledge, current gravity, 
which might indeed deviate from Einstein's one, ought to have its origin in the Planck era and 
come from a more fundamental quantum gravity/string/M-Theory approach. 

Let us briefly review the classical modified gravity theory which depends only on scalar 
curvature: 



d^xV^fiR) . (1.1) 



Here = 16ttG and f{R) is, in principle, an arbitrary function. By introducing the auxiliary 
fields A and B, one may rewrite the action ()1.1|1 as 

/ = J_ / d^^^ [B {R-A) + f{A)] . (1.2) 



Using the equation of motion and deleting B one gets to the Jordan frame action. Using the 
conformal transformation 5^,^ — > e'' g^^i, with a = — ln/'(^), we obtain the Einstein frame action 
(scalar-tensor gravity), as follows: 
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Via) = e^g (e^) - e^^/ (, (e^)) = 77^-;^- (1-4) 

Note that two such classical theories, in these frames, are mathematically equivalent. It is known 
that they are not equivalent however at the quantum level (ofF-shell), due to the use of different 
parametrizations. Even at the classical level, the physics they describe seems to be different. 
For instance, in the Einstein frame matter does not seem to freely fall along geodesies, what is 
a well established fact. 

As an interesting and specific example of the general setting above, the following action 
corresponding to gravity modified at large distances may be considered ^ 

1=^1 (fx^{R- iiR-'') . (1.5) 

K J 

Here ^ is (an extremely small) coupling constant and n is some number, assumed to be n > — 1. 
The function f{A) and the scalar field a are 



f{A) = A- iiA-"^ , cj = -ln(l + n/iA~"~^j . (1.6) 

The Friedman-Robertson- Walker (FRW) universe metric in the Einstein frame is chosen as 

3 2 

dsl = -dtl + al{tE)Y.[dx^ . (1.7) 

i=l 

If the curvature is small, the solution of equation of motion is |14j 

aE^t^ , a = — —In—, = — 7 -r 1 — • (l-»j 

n + 2 tEo 2(l + i)(n^)^ 

The FRW metric in the Jordan frame is 

3 2 

ds^ = -dt^ + a^{t) J2 {dx') , (1.9) 

i=l 

where the variables in the Einstein frame and in the physical Jordan frame are related with each 
other by 

t= e^dtE , a = e2aE, (1-10) 



which gives t ~ i]^^^ and 

(n + l)(2n + l) + 7n- 1 , , 

a~t "+2 , w = -, —, -. (l.ll) 

3(n + l)(2n + l) ^ ' 

The first important consequence of the above Eq. ()1.11|) is that there is the possibility of acceler- 
ated expansion for some choices of n (a kind of effective quintessence). In fact, if n > or 
— l<n<— ^,it follows that w < and > 0. This is the reason why such a theory ^ and 
some modifications thereof El El El E] have been indeed widely considered as candidates 
for gravitational dark energy models. The black holes and wormholes in such models have been 
also discussed il9i iSOi . 
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When —l<n<—^, one arrives at w < —1, i.e. the universe is shrinking. If we replace 
the direction of time by changing t by —t, the universe is expanding but t should be considered 
to be negative so that the scale factor a ought to be real. Then there appears a singularity at 

2 

t = 0, where the scale factor a diverges as a ~ (— . One may shift the origin of time by 
further changing —t with tg — t. Hence, in the present universe, t should be less than ts and a 
singularity is seen to appear at t = tg (for a discussion of this point, see [H]): 

ar^ {ts- t)^(^ . (1.12) 

That is the sort of Big Rip singularity. We should note that in the Einstein frame ()1.7() . the 
solution (|1.8j) gives w in the Einstein frame as 

WE = -l + ^f^>-l. (1.13) 
9{n + 1)^ 

Therefore, in the Einstein frame, there is no singularity of the Big Rip type. In general, for 
the scalar field with potential U{(p) and canonical kinetic term, the energy density p and the 
pressure p are given by 

p=^^^ + U{ip), p=^^^-Ui^). (1.14) 

Therefore w is given by 

which is bigger than —1 if U{lp) > —^Lp"^. Therefore in order that the phantom with w < —1 
is realized, one needs a non-canonical kinetic term or a negative potential. For the action (|1.3() . 
the sign of the kinetic term for the scalar field a is the same as that in the canonical action. 
Therefore, in order to obtain a phantom in the Einstein frame, at least, we need V{a) < 
or Af'{A) < f{A). For the case (|1.5j) . if p < 0, V can be negative. From (|1.6j) . however, if 
the curvature i? = ^ is small, a becomes imaginary if n > 0, what indicates that the curvature 
cannot become so small. This is quite general. Let us assume that for a not very small curvature, 
f{R) could be given by the Einstein action, f{R) ~ R; then, f'{R) ~ 1 > 0. On the other hand, 
if we also assume that when the curvature is small, f{R) behaves as f{R) ~ pR^^\ where p < 
and n > 0, then f'{R) < for small R. Hence, f'{R) should vanish for finite R, which tells 
us that a diverges since a = — In /'(A). Therefore R cannot become so small. Even in the 
case when > n > — 1, Eq. (|1.8|) tells us that and therefore a becomes imaginary. Then 
a negative p could be forbidden. At least the region where we have effective Einstein gravity 
does not seem to be continuously connected with the region where curvature is small for the 
case when a negative potential from modified gravity follows. Thus, classical modified gravity 
mainly supports the accelerated expansion. 

The paper is organized as follows. In the next section the classical dynamics of f{R) gravity in 
de Sitter space is considered. Section three is devoted to the calculation of the one-loop effective 
action in f{R) gravity in de Sitter space. The explicit off-shell and on-shell effective action is 
found. The study of quantum-corrected de Sitter geometry for specific model of modified gravity 
is done in section four. It turns out that quantum gravity corrections shift the radius of de Sitter 
space trying to destabilize it. The calculation of the entropy for de Sitter space and constant 
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curvature black holes is done in section five. Some remarks about black hole nucleation rate 
are done. The Discussion section gives some summary and outlook. In the Appendix A the 
black hole solutions with constant curvature are explicitly given for f{R) gravity. Appendix B 
is devoted to the details of the calculation of functional determinants for scalars, vectors and 
tensors in de Sitter space. 



2 Modified gravity models 

Recall that the general relativistic theory we are interested in is described by the action 

/ being a function of scalar curvature only. 

We have many possible simple choices for the Lagrangian density f{R)- Here we consider 
the following three examples: 

f (R) = R + pR^ -2A, (2.2) 
that is Einstein's gravity with quadratic corrections; 

f{R) = R-t^, ^i>0, (2.3) 

the model proposed in Ref. ^ and its trivial generalization 

/(i?) = i?-^-A.2. (2.4) 

Here we shall be interested in models which admit solutions with constant 4-dimensional cur- 
vature R = Rq, an example being the one of de Sitter. The general equations of motion for the 
model described by Eq. 1)2. 1|1 are (see, for example, |U) 

f'iR)R^u - lf{R)gf,u + (V^V, - 5m.v2) f'{R) = , (2.5) 

f'{R) being the derivative of f{R) with respect to R. If we now require the existence of solutions 
with constant scalar curvature R = Rq, we arrive at 

f'{Ro)R^, = l^g^,. (2.6) 

Taking the trace, we have the condition |21j 

2f{Ro) = Rof{R^) (2.7) 
and this means that the solutions are Einstein's spaces, namely they have to satisfy the equation 

P _ f{Ro) ^ _ Ro 

-ft^i^ — 2f'(^RQy^^ ~ 4 '^'^ ' ^ ' 

Rq being a solution of Eq. H2.7|) . This gives rise to the effective cosmological constant: 
» _ f{Ro) _ Ro 
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For the model defined by Eqs. (2.2) and ()2.3|1 we have, respectively, 

i?o = 4A, Ae// = A, (2.10) 

i?o = 3/Ui, Ae// = ±v^- (2-11) 

It is clear that this class of constant curvature solutions contains the 4-dimensional black hole 
solutions in the presence of a non vanishing cosmological constant, like the Schwarzschild-(anti)de 
Sitter solutions and all the topological solutions associated with a negative Ae//. In particular, 
with Ae// > 0, there exist the de Sitter and Nariai solutions, while for Ag// < there exists the 
anti de Sitter solution. For the sake of completeness, in Appendix^we shall review all of them. 



3 Quantum field fiuctuations around the maximally symmetric 
instantons 

In this Section we will discuss the one-loop quantization of the model on the a maximally 
symmetric space. Of course this should be considered only an effective approach (see, for instance 
[Zj). To start, we recall that the action describing a generalized Euclidean gravitational theory 
is 

^e[9] = -y^ I d^x^fiR), (3.1) 
where the generic function f{R) satisfies -on shell- the condition 

f'iRo) = ^ , (3.2) 
no 

which ensures (as we have seen) the existence of constant curvature solutions. In particular, we 
are interested in instantons with constant scalar curvature Rq, being also maximally symmetric 
spaces, namely having covariant conserved curvature tensors. An example is the 5(4) de Sitter 
instanton. For maximally symmetric space, we have 

o(0) _ ^0 / (0) (0) „(0)^(0)\ „(0) _ ^0 (0) p_ R _ 12 , . 

Another symmetric background is -ff(4), associated with a negative cosmological constant. For 
the <S'(2) X S'(2) instanton, the consideration we are going to extract are not valid, because it is 
not a maximally symmetric space and Eq. (|3.3)) . which we shall use several times from now on, 
does not hold true. In that case, one should make use of the techniques described in refs. [221 US]- 
To start with, let us consider small fluctuations around the maximally symmetric instanton 



9^J = 9if + , g'^ = - h^' + h^'K + Oih^) , h = g^^^^^h, . (3.4) 
As usual, indices are lowered and raised by the metric g^^^ . Up to second order in hij, one has 



^ -.l + \h + \h^ --hijW^ +0{h') (3.5) 



2 8 4 



5 



and 

R ~ Ro-^h + ViV.h'^ -Ah 
4 

+^ h^'^hik - T Vi - ^ Vkhi.V^h'^ + VihiVih^'' - - V^hikV'h^^ , (3.6) 

where represents the covariant derivative in the unperturbed metric gf^\ 

By performing a Taylor expansion of f{R) around Rq, again up to second order in hij, we 

get 



SE[g] 



1 



where, up to total derivatives 
1 



-9 



(0) 



fo + -{2fo-Rof^)h + C2 



(3.7) 



2 /o ^ij ^i^i^rVs /irs 



1 



h [+4i^'A - 2/^ + iio/o ] V.V,- /i. 



48/^A 2 - 24 (/^ - i?o/o ) A + 12/o - 8i2o/o + 3/?^/^ 



(3.8) 



For the sake of simplicity, we have used the notation /o = f{Ro), /q = f'iRo) and /q = f"{Ro), 
and in what follows we will set X = \{2fQ — i?o/o) (deviation from the on-shell condition). 

It is convenient to carry out the standard expansion of the tensor field hij in irreducible 
components namely 



hij = hij + Vi^j + Vj^i + ViVjCJ + -gij{h - Aqcj) , 



(3.9) 



where a is the scalar component, while and hij are the vector and tensor components with 
the properties 







hii = , 



(3.10) 



In terms of the irreducible components of the hij field, the Lagrangian density, again disregarding 
total derivatives, becomes 



jC-2 = ^h'^{3f[,A2-3fo + Rofo)kj 



+- (2/o - i?o/^) r (4Ai + i?o)?i 



32 
1 

+ 32^ 



16 



2 ^// 




9/^ A 2 - 3(/^ - 2i?oi^' ) A + 2/o - 2i?o/^ + i?^/, 
9/^A^-3(/^-2i?o/^)A3 
-(6/o - 2i?o/o - ^o/o ) A g - Ro{2fo - Rofo ) A o 



9/^ A 3 + 3(/^ - 2Rof;^ ) A 2 + i?o(/^ - i?o/o ) A n 



(3.11) 
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where A o, A i and A 2 are the Laplace-Beltrami operators acting on scalars, traceless-transverse 
vector and traceless-transverse tensor fields respectively. The latter expression is valid off-shell. 
In obtaining such expression, due to the huge number of terms appearing in the computation, 
we have used a program for symbolic tensor manipulations. 

As it is well known, invariance under diffeormorphisms renders the operator in the {h, a) 
sector not invertible. One needs a gauge fixing term and a corresponding ghost compensating 
term. We consider the class of gauge condition, parametrized by the real parameter p, 



Xk 



jrijk 



1 + 



the harmonic or de Donder one corresponding to the choice p = 1. As gauge fixing, we choose 
the quite general term [Tj 



9f 



Gij = agij + l3gijA , 



(3.12) 



where the term proportional to a is the one normally used in Einstein's gravity. The corre- 
sponding ghost Lagrangian reads [Zj 



gh 



6 

6e^ 



(3.13) 



where and Bj. are the ghost and anti-ghost vector fields, respectively, while 6 is the 
variation of the gauge condition due to an infinitesimal gauge transformation of the field. It 
reads 



6 hij = ViEj + VjEi => 

Neglecting total derivatives, one has 
Cgh = B' {a Hij + (3 A H,,) & 
where we have set 



5 y* 1 — 



gij A + 



(3.14) 
(3.15) 
(3.16) 



In irreducible components, one obtains 

2 



-9/ 



a 
2" 



4 / 



16 



hAoh 



16 



^ Ao + 



/i (Ao + ^l Ana 



Rr 



Aoa 



2 



Rq 



3p 



e ( Ai + ^i Ai^ + ^/i 



Rq 



R 



R 



+ - Ao + - Aoa 



16 



^/i Ao + -^ Ao/i-T-:a Ao + -^ Ao + ^ Aq^ 



16 



Ro 



Rq 



(3.17) 
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gh 



Rq 



C3 



Rq 



+ 



p-3 



2 

Rq 



b A 



■ 



b Ao + -^ A 



Rq 



Rq 
p-3 



Aoc 



AnC 



(3.18) 



where ghost irreducible components are defined by 
Ck = Ck + Vkc, VkC' = 0, 

Bk = Bk + Vkb, VkB^ = Q. (3.19) 

In order to compute the one-loop contributions to the effective action one has to consider the 
path integral for the bilinear part 

C=C2+ Cgf+ Cgh (3.20) 

of the total Lagrangian and take into account the Jacobian due to the change of variables with 
respect to the original ones. In this way, one gets jlUl [7| 



= {det Gij)'^/^ J D[hij]D[Ck]D[B''] exp J d'^x^C 



(det Gij )"^/^ det Jf ^ det 



1/2 



D [h] D [hi,]D [e]D [a] D[Ck]D [B^] D [c] D [b] exp 



d X ^/g C 



(3.21) 



where J\ and J2 are the Jacobians due to the change of variables in the ghost and tensor sectors 
respectively ^0]. They read 



Ji = A 



■ 



J2 



A ^0 



R 



■0 



and the determinant of the operator Gj^, acting on vectors, can be written as 

a ' 



det Gj, = const det A 1 + 



, ^ Rn a 
det( Ao + ^ + -^ 



(3.22) 



(3.23) 



while it is trivial in the case /3 = 0. 

Now, a straightforward computation leads to the following off-shell one-loop contribution to 
the "partition function" 



_r(i) 



= Z« = det 



Ai + 



X det 



/?A: 



X det 



(/5-3)Ao-i?o 



X det 



A< 



i?0 
~6~ 



X det 
X det 



2(a + /?Ai) Ai + 



X 
Wo 
Ro 



1/2 

X det 
Rq 



/?(Ao + ^)+a 

-1/2 



1/2 



A2 + 
+ X 



3 

-1/2 



/^(Ao + f 



2fo_ 
3Ro 



+XCi + X'C2 



/J^Ao + ^l+a 

-1/2 



ip-3)Ao-Ro 



(3.24) 
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where F^^^ is the one-loop contribution to the partition function and Ci and C2 are operators, 
which read 



Co 



2/fio 2aRl (3Rl fj^Rl 

3 6 3 

4/ ^ ^ n[3Rl bf'^Rl 

-— +3aiio H — - + — — - 

3 12 3 

2aRop P RqP uRqp^ 



3 6 3 

/ 15/3i?o , „ w, p „ 7PRop ap 
+ 3 a H h2/o Ro-2ap h 



6 




(Ao + i?o). 



(3.25) 



Equation (|3.24j) reduces to the corresponding one in Ref. JHl) when /" = (Einstein's gravity 
with a cosmological constant). For another approach to the same problem see Ref . |24| . 

In the derivation of (|3.24|) . it is understood that the functional determinants has been regu- 
larized by means of zeta function regularization (see, for example jl2l \1'6\ ). However, we should 
remind that within the zeta function regularization, it is no longer true that 

det = det ^ det B , (3.26) 

where A and B are two (elliptic) operators. In fact, in general, one has 

detAB = e''^^'^^detAdetB, (3.27) 

where a{A,B) is a local functional called multiplicative anomaly (see, for example, |25| I26j). 
As a consequence, in the above manipulations, we have assumed the multiplicative anomaly 
to be trivial, namely a{A, B) = 0. This is justified since here we are limiting ourselves to the 
one-loop approximation, and in such a case, a non-trivial multiplicative anomaly, which is a 
local functional of the fields, may be absorbed into the renormalization ambiguity |27j . 

Furthermore, another delicate point should be mentioned. The Euclidean gravitational ac- 
tion, due to the presence of i?, is not bounded from below, since arbitrary negative contributions 
can be induced on R, by conformal rescaling of the metric. For this reason, we have also used 
the Hawking prescription of integrating over imaginary scalar fields. Furthermore, the problem 
of the presence of additional zero modes introduced by the decomposition (|3.9|) can be treated 
making use of the method presented in Ref. jlUj . 

As one can easily verify, in the limit X — > (on-shell condition), Eq. (|3.24)) does not depend 
on the gauge parameters and reduces to 



^on-shell — Isigo) + F^,^^ 



-shell 



-A2 + 



Rr 



Indet 

2 

1 , , 

+- Indet 



Ri 



Rq 



-A, 



I 2/0 
3 3iio/o 



(3.28) 
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As usual, an arbitrary renormalization parameter £^ has been introduced for dimensional reasons. 
When /g = 0, namely in the case of Einstein's gravity with a cosmological constant, f{R) = 
R — 2A, one obtains the well known result [Hj 1^1 El 



^on-shell = -^^(fi'o) + T 



(1) 

on— shell 



1 



In det 



-A 



A ^0 

6 

4 



(3.29) 



In order to simplify the off-shell computation, we choose the gauge parameters p = l,/3 = 
and a = 00 { Landau gauge). Thus, we obtain 



24vr , 1. . 
— "2 /o + - In det 



GRi 
1 







-A. 



R^ X + 2fo 
6 A -2/0 



In det 



+-ln det 



Ai- 
-Ao 



4 



1 



In det -A 



'0 



^0 
2 



5Ro X - 2/0 



12 



5Ro , A -2/0 



12 



+ 



6 



X-fo 
3/^ 



(3.30) 



Recall now that the functional determinant of a differential operator A can be defined in terms 
of its zeta function by means of (see for example ^1 E] ) 



Res > 



D 



(3.31) 



lndet(^2^) = -C'{0\fA) = -C'{0\A) + lnfC{0\A) , (3.32) 

where the prime indicates derivation with respect to s. Looking at Eq. (|3.80j) . we see that the 
one-loop effective action can be written in terms of the derivative of zeta functions corresponding 
to Laplace-like operators acting on scalar, vector and tensor fields on a 4-dimensional de Sitter 
space. In all such cases, the eigenvalues of the Laplace operator are explicitly known and the 
zeta- functions can be computed directly using Eq. 1)3. 31(1 . For the reader's convenience, we have 
reported in the Appendix El all the details of the method used in the explicit computation for 
the example that will follow. 

Finally, equations (IOn|) . (HO^ and (|09|l . (f05|l . and (|R38|l in Appendix El lead to the 
off-shell one-loop effective action 

r= ^/o-^Q2(a2) + iQi(ai) + ^Qo(ao)-^Qo(a+)-^Qo(a-), (3.33) 
where (see App. El) 
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02 = 


-j + q2, 


g2 = 2 




13 25 




cti = 


T + ^^ = T' 


qi = 




9 33 


qo = Q 


"0 = 


4+'^° = T' 



X + 2fo 
X-2fo 



(3.34) 
(3.35) 
(3.36) 
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9 X -2fo 

a± = - + q±, q± = 5 + 2 




24(1 + 2^). (3.37) 

Now, we would like to present the explicit example for the model described by Eq. ()2.4() . 
First of all we consider the simplest case in the class of models defined by Eq. (|2.4|) . thus 

f{R) = R-^-^i2, X = flo - 3^ - 2/i2 (3.38) 

We may eliminate X and get 

57^1 +i?o (32^2 -7i?o) ,„ 
= ' ^ ^ 

33 Rn 



a± = — + ^±—jR^-36fil + l2fiiRo{Ro-2fi2)- (3.40) 
4 /ii f^i^ 

Hence, the one- loop effective action in f{R) gravity in de Sitter space is found. 

In the next section the above effective action will be applied to study the back-reaction 
of f{R) gravity to background geometry. However, several important remarks are in order. 
As usually, any perturbative calculation of the effective action in quantum gravity is gauge 
dependent. The way to resolve such a problem is well-known: to use the gauge-fixing independent 
effective action (for a review, see jTj)- More serious problem is related with the fact that 
quantum gravity under investigation is not renormalizable. Then, generally speaking, higher 
order corrections are of the same order as one-loop ones (the same is applied to all previous 
quantum considerations of Einstein gravity). As a result all one- loop conclusions are highly 
questionable as they maybe spoiled by higher loops effects. In this respect, the results of our work 
are definitely useful in the following sense. One can expect that perturbatively renormalizable 
gravity maybe constructed for some version of f{R) gravity. (So far only higher derivative 
gravity is known to be renormalizable). In this case, our work gives necessary background for 
one-loop quantization of such theory. From another side, to get the meaningful results with 
non-renormalizable quantum gravity one may apply the exact renormalization group scheme. In 
such a case, higher loop effects are not important. Our work maybe considered also as necessary 
and important step in this direction. Indeed, it is technically clear how to construct the exact 
RG equations for f{R) gravity using results of this section in the analogy with Einstein gravity 



4 Quantum-corrected de Sitter cosmology 

Let us consider the role of quantum effects to the background cosmology. So far, such study 
has been done for Einstein or higher derivatives gravity only. In order to see the difference with 
such models, we take the example of modified gravity with the action 

f{R)=R-i^. (4.1) 

It is interesting to investigate the region where curvature is not very big, as otherwise the 
classical theory is effectively reduced to Einstein's gravity. Moreover, if curvature is small one 
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can neglect the powers of curvature in the one-loop effective action supposing that logarithmic 
terms give the dominant contribution. The parameter fii is chosen to be very small in order to 
avoid conflicts with Newton's law. As a result, one obtains 

r,fl,) . ^ f iio - i^U fa + ^) m f ^ V (4.2) 



GRl V RoJ \ m+RlJ \ 12 

Here a and /3 are constants. It is assumed that the curvature is constant, R = Rq. Let us find 
the minimum of F with respect to Rq. One can write F' (i?o) as 

r'{Ro) = F{Ro)-G{Ro) , 
247r / . 3^1 \ 



When 

Ro = Rc = ^Mn , (4.4) 

F(i?,o) = 0, which corresponds to the classical solution. When Rq ~ 0, F (Rq) behaves as 
F {Ro) - and when Rq ^ +oo, F (Rq) ^ Since 

there is a minimum for F (Rq) when Rq = -v/6/ii > Re- On the other hand, if a 7^ 0, G{Ro) 
behaves as 

G{RQ)^-^(a + ^) , (4.6) 



Ro \ 111 
when Rq ^ ^ and 

G(i?o) - , (4.7) 

when i?o ^ +cxd. Hence for i?o > 0, if a < 0, r'(_Ro) > when i?o ^ +0 and T'{Rq) < when 
i?o — +00. Therefore, there is a solution which satisfies F'(i?o) = if a < 0. When a> 0, the 
existence of the solution depends on the details of the parameters. 
In case a = 0, when Rq — > 0, G{Rq) behaves as 

G{Ro) - --4r ' (4-8) 
and when Rq +00, we find 

«<^»)-|'"(^)- 

When Rq> 0,if (3 > 0, r'(i?o) > when Rq +0 and T'{Rq) < when Rq +00. Hence 
even if a = 0, when /3 > 0, there is a solution for equation r'(i?o) = 0. When /3 < 0, the 
existence of the solution depends on the details of the parameters again. Thus, there could be 
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a positive non-trivial solution for Rq, which describes the quantum-corrected de Sitter space. 
One may play with the parameters of the theory under consideration in such a way that the 
quantum-corrected de Sitter space can provide a solution to the cosmological constant problem. 
The above results indicate that the classical de Sitter solution (|4.4j) can survive when one takes 
into account the quantum corrections. A similar consideration can be done for any specific f{R) 
gravity. 

Let us demonstrate that indeed with some fine-tuning the obtained effective action maybe 
used to resolve the cosmological constant problem. One can present (|3.33|) corresponding to 
(|3U) as 



Ro-j^)+Q{l';Ro,f^i) . (4.10) 



r = ^(R,-l^ 

GRq V ^0 
In general Q (p; Rq, fJ-i) has a structure as 

Q(/^i?o,m)=Qof^V"^ + Qi( — I • (4-11) 




By the condition that T takes a minimum value with the variation over Rq, we obtain 

^ ^ ^ ^ 247r fj_ ^ dQ {l^;Ro,fsi) _ ^^^^2) 

QRq G \ Rq Rq J QRq 

The convenient choice between the parameters is 
/12\2 2 



y^^j (4.13) 
Here cq is a constant which could be determined later. Then Q has the following form: 



g = Q ^ =Qo ^ In -J^ UQi U^ . (4.14) 




We now consider the possibility that the vanishing cosmological constant could be obtained 
by (fine-) tuning the parameters. The corresponding condition that the vacuum energy, or 
cosmological constant, vanishes requires 

r = , (4.15) 

which may be solved with respect to fii as ni = /ii (Rq), what gives 

dr dfii dr 



By combining (|4.12|) and (|4.16j) with (|4.14)) . one gets 




which gives 




247r/if 
~GR^ 



(4.18) 
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Then by using KT^ . (EHH), and KWi . it follows 



Hence, 

32 



i?o = /^i • (4.20) 
By using and we find 

Q(1) = 0. (4.21) 



Then Eq. (|4.14|) shows that 

Co = e «o{i) . (4.22) 

Therefore, including the quantum corrections and (fine-)tuning the theory parameters, we may 
obtain the solution expressing the vanishing (effective) cosmological constant. Of course, such 
solution of cosmological constant problem is one- loop, and in higher orders better fine-tuning 
maybe required. 

The effective action ()4.2() has been evaluated in the Euclidean signature, in which case we 
should recall that the 4d de Sitter space with positive constant curvature Rq becomes a sphere 
of radius 

/ 1 o 

(4.23) 

The volume (area) of the sphere V is 
SiT^a* 384^2 

Identifying / d^x^ ~ y = SSAtt^RI one may reasonably assume the local effective Lagrangian 
corresponding to (|4.2j) to be 

VgLesiR) , (4.25) 




3847r2 



The effective equation of motion is 

= l9^.uL,s{R) - R^,,L'^s{R) + ^f^^uK^iR) - g,,,V^L'^s{R) , (4-26) 

with the curvature being covariantly constant, V pR^y = 0, Eq. (|4.26|) reduces to T'{Rq) = in 
(|4.3I) . Supposing the FRW metric with flat 3-dimensional part, 

d52^_^^2^g2aW ^ ^ ^^^^^^ 

1=1,2,3 
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the (t, t)-component of H4.26() has the following form 
1 







^Lcff (m + 12hA +?,{h + (qH + 12/^2 



(4.28) 



The Hubble parameter H is defined by i/ = ^, as usual. We now split Lcs{R) = Lc{R) + Lg{R), 
with 

^^^^^ = -^[iirr2+iR ' 



LgiR) 



R^ [a + 



In 



fR 
12 



(4.29) 



Let us assume that Lq[R) are much smaller than Lc{R) and consider the perturbation from the 
classical solution in (|4.4j) . by putting 



' Rc / \/3mi 



^'12 V 12 



or 



R = R^ + 5R , 5R = 65h + 24.hc5h . 
Note that Lc{R) contains the quantum correction. From (|4.28|) it follows that 

= -18/icL"o(5/i - 54hlL%6h + (-6/1^^0 + 72/icL"o) 6h - ^Lgo + Sh^Lg 



(4.30) 



(4.31) 



qO 



247r •• 72tt ■ 



2887r/ic 
G 



6h 



+2Ak 
2Att 



a 



768/xi 



+ I2hl ( a + ^ 



Gh 



5h + 3hc6h - 12hi6h 



2 a 



In ( Z3/i2 ) + ( a + 



192//1 

P 



Here L' 



with 



cO 

6h -- 



hi 



27r I V 768/xiy V ^7 V 192/ii 
= L'^ (iio) and L"q = (i?o)- Then the solution is 

/ii +^+e°+* + ^_e^-* , 



(4.32) 
(4.33) 



A, 



Ghc 
247r 
1 
2 



2 a 



768;Ui 
3 ± \/57l /ic . 



In ( ) + ( a + 



P 



192/ii 



(4.34) 



Since > 0, the de Sitter solution becomes unstable under the perturbation. 

Thus, for a specific modified gravity model we have demonstrated that the quantum gravity 
correction shifts the radius of the de Sitter space trying to destabilize the de Sitter phase. This 
may find interesting applications in the study of the issue of the exit from inflation or in the 
study of the decay of the dark energy phase. 
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5 Black hole nucleation rate 



We have remarked (see Appendix ^ that within the modified gravitational models we are 
dealing with, there is room for black hole solutions, formally equivalent to black hole solutions 
of the Einstein theory with a non vanishing cosmological constant. As in the Einstein case, 
one is confronted with the black hole nucleation problem We review here the discussion 
reported in refs. [221 123]. 

To begin with, we recall that we shall deal with a tunneling process in quantum gravity. 
On general backgrounds, this process is mediated by the associated gravitational instantons, 
namely stationary solutions of Euclidean gravitational action, which dominate the path integral 
of Euclidean quantum gravity. It is a well known fact that as soon as an imaginary part appears 
in the one-loop partition function, one has a metastable thermal state and thus a non vanishing 
decay rate. Tipically, this imaginary part comes from the existence of a negative mode in the 
one-loop functional determinant. Here, the semiclassical and one-loop approximations are the 
only techniques at disposal, even though one should bear in mind their limitations as well as 
their merits. 

Let us consider a general model described by f{R) with Ag// > 0. Thus, we may have de 
Sitter and Nariai Euclidean instantons. Making use of the instanton approach, we have for the 
Euclidean partition function 

Z ~ Z{S4) + Z{S2 X S2) = (54)6-^(^4) + Z(^)(52 X 82)6-^^^^''^'^ , (5.1) 

where / is the classical action and Z^^^ the quantum correction, typically a ratio of functional 
determinants. The classical action can be easily evaluated and reads 

I{S,) = -P^, I{S2xS2) = -^. (5.2) 

At this point, we make a brief disgression regarding the entropy of the above black hole 
solutions. To this aim, we follow the arguments reported in Ref. [TP. If one make use of the 
Noether charge method for evaluating the entropy associated with black hole solutions with 
constant curvature in modified gravity models, one has 

S = ^fiRo) ■ (5.3) 

As a consequence, in general, one obtains a modification of the "Area Law" . For stable models 
like (|2.2|) with p > (see below), one has /q = 1 -|- 8pA. For the model p.3|) . /q = |, and thus 

m. 

In the above equations. Ah = 47rr|^, being the radius of the event horizon or cosmological 
horizon related to a black hole solution. This turns out to be model dependent. It is interesting 
to note that for unstable modified gravity (with negative first derivative of /) the entropy may 
be negative! 

We are interested in the case of de Sitter space, and we have 

127r 487r , , 
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Thus, for the de Sitter solution, 
1 Itt 

S{S,) = ^f'{Ro) (5.6) 

and if we take into account the Eq. ()2.7() , one gets 

SiS^) = -I{S^) . (5.7) 

With the help of the one-loop effective action one can calculate quantum correction for classical 
entropy. 



We may introduce the free energy F = where /3 = 27ry is the inverse of the Hawking 
temperature for the de Sitter space. As a consequence, we have [23] 

Z~Z«(54)e-^-^, (5.8) 

where 

^ = ^^^^^ - ms.) ■ ^'-'^ 

The rate of quasiclassical decay in the de Sitter space is present as soon as F has a non vanishing 
imaginary part and it is given by = 2 Im F . When /(i?) = R — 2A, the Einstein case, it turns 
out that Z({S2 X 5*2) has an imaginary part but Z{S4) is real. As a result, in the Einstein case, 
the nucleation rate is j22| \2'A\ 

Within our generalized models, the dynamics of the gravitational field is different. In fact, also 
in the de Sitter case, due to the presence of an additional term in the on-shell one-loop effective 
action related to the operator Lq, see Eq. (|;-{.28j) . there exists the possibility of negative modes. 
In fact, from Appendix B, one has 

12 r + ^Rof!^ 

where n = 0, 1, 2, 3, ... It is clear that we have negative modes as soon as 

4^<^- (5-12) 

For example, for the model 

f{R)=R-^, (5.13) 

the quantity 1)5. 12(1 is always negative and one obtains, at least, two negatives modes. 
For the model 

f (R) = R + pR^ - 2A, (5.14) 

there are no negative modes as soon as p > 0, in agreement with the classical stability observed 
in ref. |2^. For p < 0, one has only a negative mode when 



A,.(io) = i^('!' + 3n-4 + -ia. , (5,11) 



P<-^. (5.15) 

Thus, in general, within these specific modified gravitational models, de Sitter space is unstable 
due to quantum corrections. 
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6 Discussion 



In summary, we have here calculated the one-loop effective action for general f{R) gravity in 
the de Sitter space. Generalized zeta regularization has been used to obtain a finite answer for 
the functional determinants in the effective action, what has proven to be a very convenient 
procedure. The important lesson to be drawn from this calculation, generalizing the previous 
program for one-loop Einstein gravity in the de Sitter background, is that quantum corrections 
tend to destabilize the classical de Sitter universe. The constant curvature black hole nucle- 
ation rate can be discussed within our scheme. Note that, typically, f{R) models may contain 
(depending of course on the explicit form of /) more negative modes, aside from the Einstein 
one. 

Another lesson we have learned is that the inverse powers of the curvature (which are im- 
portant at current epoch) do not arise from the perturbative quantum corrections, as has been 
explicitly demonstrated here. Some remarks about black holes and associated entropy in modi- 
fied gravity have been also made. 

There are several directions in which our results can be extended and applied. First of all, 
one can repeat the whole calculation of the functional determinants for the case of Anti-de Sitter 
space. Then, the one-loop effective action found in the third section can be given also for the 
Anti-de Sitter universe. In relation with the AdS/CFT correspondence, our approach can be very 
interesting for the study of the (de) stabilization of such a universe against quantum corrections. 
On the other hand, it can be also important in the search for the solution of the cosmological 
constant problem within hyperbolic backgrounds in the large distance limit 30 , what will be 
discussed elsewhere. Second, having in mind the current interest on gauged supergravities in de 
Sitter space, one can try to formulate the theory of f{R) gauged supergravity and to study the 
relevance of quantum effects for such theories. 
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A Black hole solutions in modified gravity 



Here we revisit the class of general static neutral black hole solutions in four dimensions and 
non vanishing cosmological constant. Let us start with the usual Ansatz for the metric 



ds'^ = A{r)dt^ 



k ' 



(A.l) 



where A; = 0, ±1 and the horizon manifolds are Si = S^, the two dimensional sphere, Sq = T'^, 
the two dimensional torus, and S_i = H'^/T, the two dimensional compact Riemann surface. 
The scalar curvature and the non vanishing components of the Ricci tensor read 



R 



1 



r'^A" + ArA' + 2A-2k 



(A.2) 
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R,^ = -l-[rA" + 2A',] (A.3) 

Rrr = ^[rA" + 2A' ,] (A.4) 

Rab = gab [k -rA'-A .] (A.5) 
If we look for a constant curvature solution, we should have 

r'^A" + 4rA' + 2A - 2k = -r'^Ro . (A.6) 
The general exterior solution depends on two integration constants, c and b, and reads 

Air) = \ + k-'--^r\ (A.7) 
However, this solution has to satisfy the equation of motion 

R^.u = ^9f.. . (A.8) 

As a result, it is easy to show that these equations are satisfied with 6 = and c arbitrary, and 
we have (see for instance j2HlEl) 

A{r) = k-^- ■ (A.9) 

Since A(r) > 0, we can have k = 1 only for positive Rq, and this is the Schwarzschild-de Sitter 
solution. For Rq < 0, we may have k = 1, namely the Schwarzschild-AdS solution. We may also 
have A; = 0, with a torus topology for the horizon manifold, and k = —1, with an hyperbolic 
topology for the horizon topology, the so called topological black holes [^S]. The constant c is 
related to the mass of the black hole. The de Sitter solution is obtained when Rq > and with 
c = and k = 1, the AdS solution is obtained when Rq < and with c = and k = 1. For c 
non vanishing, one has black hole solutions. These black hole solutions may have extremal cases 
and extremal limits. The extremal case exists for k = —1. For k = 1, Rq > 0, one has only 
the extremal limit of the Schwarzschild-de Sitter solution (see, for example [22j and references 
therein), and the metric reduces to 

ds"^ = (dSi + dSi) . (A. 10) 

Rq ^ ' 

This is a space with constant curvature Rq solution of Eqs. HA.8|1 . 



B Evaluation of functional determinants 

Here we shall make use of zeta function regularization for the evaluation of the functional 
determinants appearing in the one-loop effective action, Eq. 1)3. 30(1 computed in the previous 
section. 

First, we outline the standard technique, based on binomial expansion, which relates the 
zeta- functions corresponding to the operators A, with eigenvalues A„ > and A = -^(^ — 
a), with eigenvalues \n = ■§-(An — a), Q being a real constant. With this choice, A„ and a 
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are dimensionless. We assume to be dealing with a second-order differential operator on a D 
dimensional compact manifold. Then, by definition, for Res > D/2 one has 



n 

Us) ^ Cis\A) = J2X-^ = 



12 



(B.l) 
(B.2) 



where, as usual, zero eigenvalues have to be excluded in the sum. In order to use the binomial 
expansion in ()B.2|) . we have to treat separately the several terms satisfying the condition |A„| > 
\a\. So, we write 



Cais) 



Rq 
12 



k=0 



k\T{s) 



where we have set 
FJs) = 



Gis) 



E (An-a)-% 

A„<|a|; Xnf^a 

Y: x-' = as)-F(s) 

A„>|a| 



Fis)= E 

An<|a| 

F„(0)-F(0) = iVo 



(B.3) 

(B.4) 
(B.5) 



Nq being the number of zero-modes. It has to be noted that ()B.3|) is valid also in the presence 
of zero-modes or negative eigenvalues for the operator A. In many interesting cases, Fa{s) and 
F{s) are vanishing and thus G{s) = C(s)- 

As is well known, the zeta function has simple poles on the real axis for s < D/2 but it is 
regular at the origin. Of course, the same analytic structure is also valid for the function G{s). 
One has 



n=0 



s + {n-D)/2 



(B.6) 



J(s) being an analytic function and Kn the heat-kernel coefficients depending on geometrical 
invariants. In the physical applications we have to consider, we have to deal with the zeta 
function and its derivative at zero, thus it is convenient to consider the Laurent expansion 
around s = of the functions 



T{s + k)(:{s + k) 

r{s + k)G{s + k) 



h a/c + 0(s) , 

s 

— + Ofc + 0{s) , 
s 



(B.7) 
(B.8) 



bo 


= bk 


- F{0) , 


ao = ao + 7F(0) , 




h 


= bk 


= Kj:,_2k , 


ak = a-k- T{k)F{k) , 




bk 


= bk 


= 0, 


^(A:) = C(A:) - F{k) , 


k > 



l<k<^ 

- - 2 

D 



(B.9) 
(B.IO) 

(B.ll) 
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Now, from previous considerations, one obtains 

Ca{s) = 



12 



E 

P<k<D/2 



k\ 



kl 



+Fa{s) + S ^ 
k>D/2 



k 



and finally 

Ca(0) 

c;(o) 



0<k<D/2 



0<k<D/2 



k\ 



k>D/2 



k 



(B.12) 

(B.13) 
(B.14) 



7 being the Euler-Mascheroni constant. If there are negative eigenvalues then -^^(0) has an 
imaginary part, which reflects instability of the model. 

In the paper we have to deal with Laplace-like operators acting on scalar and constrained vec- 
tor and tensor fields in a 4-dimcnsional de Sitter space 50(4). In all such cases, the eigenvalues 
A„ and relative degeneracies gn can be written in the form 



Rq 
12 



An - a 



Qn = Ci {n + V')+ C3 (n -I- vY 



An = (n + l^f , 



(B.15) 



where n = 0, 1, 2... and ci, C2, ly, a depend on the operator one is dealing with. In our cases we 
have 



Rq 

12 



3 

2 ' 



Li 



Lo 



Rq 
12 



-A- 



i?0 
12 



ci 
J/ = 

ci 
v - 

Cl 



■ 12 5 C3 



5 

2 ' 



C3 = l. 



4 ' 

7 17 

2 ' 



a 



12 ' ''3 



5 

3 ' 



(B.16) 
(B.17) 
(B.18) 



where q are dimensionless parameters depending on the specific choice of f{R)- 
We note that ({s) is related to well known Hurwitz functions Ch{s, i') by 



E 9n'^n' = E [^1 + ^)^'~^ + C3 (n + 



2s-3 



n=0 



n=0 



and 



ciCh {2s - + csCh {2s - 

ciQh {2s csCh {2s -3,u)- F{s) 

ciCh {2s -l,u + n) + C3Ch {2s -S,u + h) , 



(B.19) 



(B.20) 



21 



h being the number of terms not satisfying the condition A„ > \a\. In order to proceed, we have 
to compute the quantities bk and Ofc for A: = 0, 1, 2. To this aim, we note that Hurwitz functions 
have only a simple pole at 1 and, more precisely. 



Ch{s + 1,u) = --iP{u) + 0{s), 
s 



(B.21) 



tp{s) being the logarithmic derivative of Euler's gamma function. After a straightforward com- 
putation, we get 



^0 = C(o) = ci C// (-1,1^) + csCh (-3,z^) , ^i^Y' ^2 = Y, 



(B.22) 



ao 

0-2 



C'(0)-7C(0) 

ci [2C'h (-1, u) - -fCn (-1, 1^)] + C3 [2C'h (-3, i/) - tCh (-3, i^) 

7-1" 



-ci 



ci Cj? (3, jy) - ca 



Using (|B.14|) we obtain 



k=0 
2 



kl 



In- 



12 



+ ^^j + ^o.{0) + 2^ 



fc=0 



fc=3 



k 



(B.23) 
(B.24) 

(B.25) 



(B.26) 



Now we have to consider separately the operators Lq,Li,L2 we are dealing with and explicitly 
compute bk, ak, and G{k) using (RT^I) . ((b22|)-(|E211, and (jHlBll - fRlSl) . 



The scalar case 

The eigenvalues of Lq are of the form 



\ _ ^0 
-^"-12 



3 

n -\ — 
2 



a = ^ + 



n = 0,1,2. 



(B.27) 



This case is model dependent since the parameter q explicitly depends on the choice of the 
Lagrangian f{R)- Then one could have zero modes and also negative eigenvalues, but we take 
them into account by the functions and F, both of which will in general appear in the final 
result. For k > 3, we have 



G{k) = -^CHl2k-l,^ + n] +^CH(2k-3,^ + h] , 



(B.28) 



where A„ > |a| per n > h. Then 



Qo{a) = C'aW^Lo 



17 a o? 

No \ 

^ 2880 24 12 



In- 



12 
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+ 3 [3^a(0) + 4C'^(-3, 3/2) - C^(-l, 3/2)] 



72F(1) + 11 - 6'(/'(3/2) 



a 
72 



12 F(2) + 4^(3/2) + 7Ci?,(3) - 10 



a 
24 



■ 

fc=3 



k 



(B.29) 



C_r(s) being the Riemann zeta function. One of the three scalar Laplacian-hke operators ap- 
pearing in the one-loop effective action ()3.33|) does not depend on the model since for such case 
a = ao = 33/4. Then Aq and Ai are smaller than q (n = 2) and so 



FUs) = {-6)-' + 5i-2y 



Fis) 



+ 5 



25 



(B.30) 



G{k)-- 

From these equations it follows 



5^C«(2*-4) + ^C«(2^--34 



479 £^ En 
Qo(33/4) ~ -18.32 - Gvr i + In ^ . 

We see that there is an imaginary part since there are negative eigenvalues. 



(B.31) 



(B.32) 



The vector case 

The eigenvalues of Li are of the form 



Ar 



-^0 
12 



n + 



a 



13 

a = — + 
4 



n = 0,1,2.. 



(B.33) 



For the vector case, g = 3 is a pure number and so a = 25/4 = Aq- Thus there is a zero- mode 
with multiplicity equal to 10 (A'^o = 10) and this has to be excluded in the evaluation of zeta 
function. As a consequence, we have Fa{s) = 0, F{s) = 10 a~'^, bo = bo — 10, ao = ao + IO7, 
ai = ai — 10a~^, a2 = 02 — 10 a~^, and for A; > 3 



Gik) 
Finally, 



(B.34) 



191 



Qi(25/4) ^C(0|^%) = -— In— ^ + 



''Ro 22215 



30 



12 



64 



+ 4C^(-3,5/2) 



-9C^(-l,5/2)-^C«(3)-^^(5/2) 



-18.91 



a''G{k) 

k=3 

191 



k 



In 

30 12 







(B.35) 
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The tensor case 



The eigenvalues of L2 are of the form 



^""12 



7\2 
n + - I -a 
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n = 0,1,2... 



(B.36) 



As for the scalar case, here also zero- modes could appear and/or negative eigenvalues, depending 
on the parameter q. Then, in general, we have to introduce the functions Fa{s) and -F(s). For 
A; > 3, we have 



G{k) = Cik) = Ch (2k -1,1\+^Ch hk -3 J- 



and 



Q2(a) = C(0K'L2 



8383 125a 5a'^ 



1 



In- 



F{k), 



'Rq 



(B.37) 



12 



+ 3 [3F;(0) + 20C^(-3,7/2) - 125C^,(-l,7/2)] 



72F(1) + 535 - 750^^(7/2) 



Q 

72 



324F(2) + 540V'(7/2) + 23625Cij(3) - 28486 

fc=3 



a 

648 



(B.38) 
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